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Abstract 

An existence and stability result for a class of purely radiative vacuum spacetimes arising 
from hyperboloidal data is given. This result generalises semiglobal existence results for 
Minkowski-like spacetimes to the case where the reference solution contains gravitational 
radiation. The analysis makes use of the extended conformal field equations and a gauge 
based on conformal geodesies so that the location and structure of the conformal boundary 
13^ of the perturbed solutions is known a priori. 

U 

1 Introduction 

In this article we analyse the stability of the hyperboloidal initial value problem for a class of so- 
lutions to the vacuum Einstein field equations — the so-called purely radiative vacuum spacetimes. 
j> — Purely radiative vacuum spacetimes describe gravitational radiation which comes from infinity, 

interacts with itself in a non-linear way and then disperses to infinity. In order to encode that the 
spacetime is made only of incoming gravitational radiation, one assumes that it admits a Penrose 
conformal extension such that the null generators of future and past null infinity (J? + and J?~) 
t-H are complete and that future and past timclikc infinity are represented in the conformal extension 

by two points i + and i . The points i are required to be regular — that is, the conformally 
rescaled unphysical spacetime admits a smooth extension which is regular at . More precisely, 
given a vacuum spacetime {M,g^ v ) it will be assumed that there exists a conformally related 
spacetime (M.,g^ v ) such that 

9iiv = Q 2 g^, (l) 

where is a suitable conformal factor. Null infinity, J^ ± , is characterised as the set of points for 
which 

9 = 0, d6 ^ 0. (2) 
On the other hand, the points in a purely radiative spacetime are such that 

O = 0, dO = 0, Hessian is non-degenerate. (3) 

It is perhaps worth noting that the corresponding points i in a spacetime with a black hole (say, 
the Schwarzschild spacetime) do not satisfy these conditions. 

It has been know for a long time that the question of the existence of non-trivial examples of 
purely radiative solutions to the Einstein field equations hinges very delicately on the structure 
of the spacetimes in a neighbourhood of the so-called spatial infinity, i° — see e.g. [TUJ, [T2l 
1141 115] - A way of getting around the so-called i° -problem is to consider initial value problems 
posed on hyperboloidal hyper surf aces, that is, spacelike hypersurfaces which can be thought of as 
intersecting null infinity in a transversal manner. If the initial data prescribed on the hyperboloid 



- 1—1 

X 



*E-mail address: cluebbeOqmul.ac.uk 

^E-mail address: j . a. valiente-kroonOqmul . ac .uk 



1 



is suitably close to a canonical hyperboloidal data implied by the Minkowski spacetime — that 
is, the data is a perturbation of Minkowski data — then it has been proved that the resulting 
development has a conformal structure similar to that of Minkowski spacetime [9] . More precisely, 
the future development of initial data admits a conformal compactification such that the locus 
of points, ,y + satisfying the condition ^ has the topology S 2 x [0, 1), and that there is a single 
point, i + , satisfying the conditions ([3| where the null generators of null infinity converge. The 
original proof of this result was carried out using the so-called conformal Einstein field equations 
discussed in [5J El [5] . These equations have the property of being regular even at the points 
where the conformal factor vanishes and are such that a solution to them implies a solution 
to the vacuum Einstein field equations. A property of the evolution systems implied by these 
conformal field equations is that it includes the conformal factor as one of the unknowns, and 
thus any discussion of the structure of the conformal boundary has to be performed a posteriori. 
The results of [3] have been generalised to the case of the Einstein field equations are coupled to 
Maxwell and Yang-Mills fields [TT] . 

A more general system of conformal field equations has been introduced in [13] . This extended 
conformal Einstein field equations are expressed in terms of so-called Weyl connections — torsion 
free connections (not necessarily Levi-Civita) which preserve the conformal structure. Thus, 
the field equations contain more gauge freedom. In particular, using the extended conformal 
equations it is possible to introduce gauges — conformal Gaussian gauges — which are based on 
certain type of conformal invariants the so-called conformal geodesies. The advantages of these 
conformal Gaussian systems are twofold: combined with the extended conformal field equations, 
they imply simpler systems of evolution equations; and, if working in vacuum, they provide an a 
priori canonical conformal factor so that the location and properties of the conformal boundary 
are known and can be adjusted according to need before the propagation equations are solved. 
In |27j . the extended conformal Einstein field equations and conformal Gaussian gauge systems 
have been used to provide a new version of the proof of the stability results for Minkowski and 
de Sitter spacetime. In particular, the analysis of the properties of the conformal boundary of 
the spacetimes obtained is much simpler and transparent. 

The proofs of the stability of perturbations of the Minkowski spacetime in [5] and in [57] 
make use of an adaptation of a very general result on the existence and stability of symmetric 
hyperbolic systems by Kato [2H E3 [26]. The possibility of using these theorems to directly 
prove the global existence of solutions to the Einstein field equations stems from the use of a 
conformally compactified picture where timelike infinity is at a finite position. As a result, Kato's 
theorems allow to provide very compact global existence results which require almost no PDE 
(partial differential equations) analysis. The latter should be contrasted, for example, with the 
analysis of the non-linear stability of the Minkowski spacetime carried out by Christodoulou & 
Klainerman [3] where precise asymptotic estimates of the solution are obtained. 

In the present article we generalise the methods and result of [57] to reference solutions which 
are purely radiative spacetimes in the sense discussed in the first paragraphs of this introduc- 
tion. A procedure to construct an infinite number of these purely radiative spacetimes has been 
given in [10] . The possibility of generalising the existence and stability results of the Minkowski 
spacetime using this class of solutions has been suggested in [TT] . This construction depends 
upon a remarkable connection between stationary solutions to the Einstein field equations and 
spacetimes with vanishing mass. The reference spacetimes thus obtained are only known in an 
abstract sense and without much information about their geometry at hand. It is remarkable that 
although one has very little explicit information about the reference solutions, it is still possible 
to use the methods of |27j to obtain a stability result. Our main result can be summarised in the 
following form. 

Theorem. Given a sequence of multipole moments for a static solution to the Einstein vacuum 
equations subject to a suitable convergence condition, there exists a hyperboloidal initial data set 
whose future development is a purely radiative spacetime admitting a conformal compactification 
with a future null infinity satisfying conditions ^ and a point representing timelike infinity 
satisfying conditions t3y. A small enough perturbation of the hyperboloidal initial data has a 
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future development whose conformal boundary has the same properties as that of the unperturbed 
solution. 

In order to prove this theorem one starts by considering the asymptotic end of an initial data 
set for a vacuum static solution to the Einstein field equations. This asymptotic end can be 
compactified to a 3-dimensional manifold, S, such that the point at infinity corresponds to a 
certain point i. The subsequent analysis will be restricted to a suitable neighbourhood of this 
point. The observations in |10j allow to construct an initial data set for the conformal Einstein 
field equations which is analytical in a neighbourhood of i. This initial data has the peculiarity of 
having a vanishing ADM mass. As our analysis is local to % the later feature causes no problem. 
A local existence for this data follows directly. Our interest is then concentrated to the region 
in the development which is included in (the causal past of i). As a consequence of our 

conformal setting, this region implies the existence a portion of a purely radiative spacctimc. The 
latter can be regarded as the future development of some data for a purely radiative spacetime 
prescribed on a fiduciary hyperboloid Ti.. The point i in 5 corresponds to the point i + of the 
radiative spacetime. In order to prove the existence of the reference radiative spacetime one 
employs a conformal Gaussian system based on a congruence of conformal geodesies leaving the 
hypersurface S orthogonally. 

It should be noted that this congruence will, in general, intersect the fiduciary hyperboloid in 
a non-orthogonal manner. This feature hinders a direct application of Kato's theorems in order to 
conclude the stability of the hyperboloidal data. The reason for this is that in order to make use of 
Kato's theorems one has to make use of a Gaussian system in which the conformal geodesies leave 
the initial hyperboloid with the same angle and orientation as the congruence in the reference 
solution meets the hyperboloid Ti.. It turns out that the data for such an evolution are not really 
intrinsic to the initial hyperboloid but contain information along the direction normal to the initial 
hypersurface. Thus, one needs to show, first, local existence for the perturbed hyperboloidal data 
— using, for example, a congruence of conformal geodesies which is orthogonal to the hyperboloid. 
Once this local existence has been established, one can construct initial data for the evolution 
along the tilted congruence — using a Lorentz transformation. From here Kato's theorems and 
the requirement of the perturbed data being close enough to the reference solution yields the 
required stability result. 

Outline of the article 

We start by setting our conventions in section [2] and recalling the physical and the conformal 
vacuum constraints. As we are working on S 3 and subsets thereof, we fix coordinates and vector 
fields on S 3 and briefly review the norm and the extension operator for Sobolev spaces associated 
to S 3 . In section [3] we discuss the conformal tools, which we will be using later on. These include 
conformal geodesic, their associated conformal Gaussian gauge and the conformal Einstein field 
equations. We follow the standard approach by Friedrich and derive an evolution system form 
the field equation by using space spinors. Section[4]is concerned with the construction of a purely 
radiative reference spacetime from time symmetric static vacuum data with non-vanishing mass 
given on a Cauchy surface passing through the point i at spacelike infinity. In section |4.3.3| the 
resulting radiative spacetime is analysed from the viewpoint of hyperboloidal data on a surface 
not passing through i. We note at this stage that the initial data on such a hypersurface is 
not proper hyperboloidal data due to the fact that the conformal geodesic congruence does not 
intersect the hyperboloid orthogonally. Thus we analyse the relationship between this so named 
tilted hyperboloidal data and the proper hyperboloidal data. Section [5] is concerned with the 
existence and stability result for radiative spacetimes. The structure of the conformal boundary 
is analysed prior to the evolution of the data. In section [53] we use a modified version of a theorem 
by T.Kato employed in [3] and [27] to establish the existence and stability for purely radiative 
spacetimes in theorems [3] and [l] 
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2 Basics and conventions 



We shall consider spacetimes {M.,^^) satisfying the vacuum Einstein field equations 

Rp, = o, (4) 

where denotes the Ricci tensor of a Lorentzian metric <? M „, (ji,v = 0, 1,2,3) with signature 
(+,—,—,—). The discussion of the solutions to equation ([4J will be carried out in terms of a 
conformally rescaled metric g^ v related to g^ u according to equation ([I]). In what follows, let 
V M denote the Levi-Civita covariant derivative of while R^xp, C^ u \ p , R^, R denote the 
associated Ricmann, Weyl and Ricci tensors and the Ricci scalar of <? M „. 

In the stability analysis we will compare a spacetime (_M,<? M „), evolved from some chosen 
initial data, with a given reference spacetime. Throughout the article, the latter will be denote 
by g^). We will also adapt this notation for all quantities related to the reference spacetime. 

We will make use of two types of hyperboloidal data. Proper hyperboloidal data is given with 
respect to a frame consisting of three tangent vectors to the hypersurface and its normal, while 
tilted hyperboloidal data depends on a frame not intrinsic in the hypersurface. In our setup this 
frame will be oriented along the congruence. We will distinguish between the two frames and the 
corresponding frame components, by using an underscore for proper hyperboloidal data. 

The conformal factor on the spacetime will be denoted by G, while f2 is the conformal factor 
on an initial hypersurface. 



2.1 The constraint equations 

Let S denote a timelike hypersurface of the spacetime {M,g^ v ). Let h a p, K a p denote, respec- 
tively, the first and second fundamental forms induced by the metric g^ v on S {a, 13 = 1,2,3). 
The vacuum Einstein field equations W) imply the following constraint equations on S 



f-K 2 + k ali K aS3 = 0, (5a) 
D a K aP - DpK = 0, (5b) 

where D denotes the Levi-Civita and f the Ricci scalar of the metric h a p. The first and second 
fundamental forms determined by the metrics g^ v and g^ v on S are related by 



h a p = Q?h a p, K a p — Vl (^K a p + TJi a f_ 



where Q, = ®\§ and the function £ on S denotes the derivative of & in the direction of the future 
directed g-normal of S. The traces K = h al3 K aB and K = h aP K afi are related by 

Q.K = K + 3E. 



The Hamiltonian and momentum constraints for vacuum, equations (5a) and (5b I, written in 
terms of the conformal fields f2, S, h a p and K a p read 

2flD a D a n - 3D a flD a tt + ^fl 2 r - 3D 2 - ^fl 2 (K 2 - K afj K a/3 ) + 2ttY,K = 0, (6a) 
fl 3 D a (n- 2 K aP ) - (D P K - 2&- 1 D[ j Y) = 0, (6b) 

where D denotes the Levi-Civita connection and r the Ricci scalar of the metric h a p. In particular 
if O = one has that 

S 2 + D a nD a n = 0. (7) 



In the sequel we will consider two different classes of solutions to equations ( 6a ) and ( 6b I 



The first class will consist of initial data sets for the Einstein vacuum field equations which are 
asymptotically Euclidean, and thus, the relevant initial hypersurfaces are Cauchy hypersurfaces of 
an asymptotically flat spacetime. Further details will be discussed in section [4~T| The second class 
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of solutions to be considered are hyperboloidal data sets. In this case the initial hypersurfaces 
are not Cauchy hypersurfaces of a globally hyperbolic spacetimes. Further details will be given in 
section 5.1 The character of the solutions to the conformal constraint equations (asymptotically 



Euclidean or hyperboloidal) is determined through the boundary conditions which supplement 



equations (6a I and (6b) 



2.2 Coordinates on § 3 and submanifolds thereof 

In the sequel we shall work with conformal spacetimes with time slices diffeomorphic to simply 
connected subsets of the 3-sphere, S 3 . Therefore, it will be important to have a good coordinati- 
sation of § 3 and to have a frame which is globally defined over this manifold. Following the work 
in [TO] [57] we regard S 3 as a submanifold of K 4 : 



§ 3 = i x A e 



(zl)2+(z 2 )2 + (z 3 ) 2 +(z 4 ) 4 =l} 



The restrictions of the functions x A , A= 1, 2, 3, 4 on R 4 to S 3 will again be denoted by x A . If 
we view the point N = (0, 0, 0, 1) as the North pole and V = {x 4 > 0} n S 3 as the open upper 
half sphere, then for any subset U of V with N €lA the functions (x 1 , x 2 , x 3 ) provide coordinates 
on U centered at N. The vector fields 

c 1 =x 1 d 4 -x i d 1 +x 2 d 3 -x 3 d2, (8a) 
c 2 = x 1 d 3 -x 3 d 1 + x 4 d 2 -x 2 d 4 , (8b) 
c 3 = x 1 d 2 -x 2 di+x 3 d 4 -x 4 d 3 . (8c) 

on K are tangent to § 3 and provide a globally defined frame field on S 3 . 

Given a simply connected 3-manifold TL with boundary dTL, there exists a diffeomorphism 
<f> : TL — ► U = {x 4 > 1/2} C T>, such that <j> : &H — > {a; 4 = 1/2}. Since tj> is a diffeomorphism 
we can map any tensor from TL to U and vice versa. For convenience we denote tensors on TL 
and IA by the same symbols and the underlying manifold is to be understood from the context. 
In particular we can pull-back the functions x A to provide coordinates on TL and use the vectors 



ci, C2, c 3 , as given by (8al-(8c), as a globally defined frame field on TL. 



2.3 Extending fields to § 3 

In the sequel it will be necessary to extend certain (vector valued) functions defined on a subset 
U C S 3 to the whole of S 3 . For this we follow the approach pursued in [9J. On the spaces 
C°°(U,R N ) and C 00 ^ 3 ^) of smooth M^-valued functions on U and S 3 define for m G N the 
Sobolev-like norm 




1/2 

hiu,q= i >: / i >: ia^-a^i 2 )^) , (9) 
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with Q = U, § 3 . Let H m (U,R N ) and H m (S 3 ,R N ) be the Hilbert spaces obtained, respectively, 
as the completion of the spaces C°°(U,R N ) and C°°(S 3 ,K Ar ) on the norms given by formula (|9|. 
Given w € H m (U,R N ), there exists a linear extension operator E such that 

(Ew)(x) = w(x) almost everywhere in hi (10a) 
\\Ew\\ m p<K\\w\\ mJ4 , (10b) 

with K a constant which is universal for fixed m. It is important to notice that if, for example, 
w = (hij , Xij is a solution to the conformal constraints, equations ( 6a I and ( 6b I , onW then 

j3 



Ew will not satisfy the constraints in S \U. 

3 The conformal Einstein field equations and conformal 
Gaussian gauges 

When discussing the existence of solutions to the Einstein field equations in terms of conformally 
rescaled spacetimes, it is convenient to introduce gauges based on conformally invariant structures. 
The particular conformal structures to be used in our analysis are conformal geodesies. These 
curves are autoparallel with respect to Weyl connections — that is, torsion-free connections which 
preserve the conformal structure, but are not necessarily the Levi-Civita connection of a metric. 
Conformal geodesies were first applied in the context of General Relativity in [20]. The idea of 
conformal Gaussian coordinates based on congruences of conformal geodesies has been introduced 
in [T3J [T3] . It has been shown in [T7] that conformal Gaussian coordinates can be used to cover 
even the strong field regions of the Schwarzschild spacetime. 

In this section we present a brief discussion of the properties of gauges based on conformal 
Gaussian coordinates, and of the associated propagation equations implied by the extended con- 
formal Einstein field equations. For further details we refer the reader to [T7l[27J[3T] and to the 
references mentioned in the previous paragraph. 

3.1 Conformal geodesies 

Given a spacetime (A4,g Ml/ ), a conformal geodesic x^(t) and its associated 1-form 6 m (t) are 
solutions to the system of equations 

i^V^ + 2(b, x)x fi ~ g(x, x)¥" = 0, (11a) 

x v V u b„ - (b, x)b„ + iff»(6, b)x„ = Z AA1 i A , (lib) 



where denotes the Schouten tensor of the Levi-Civita connection V, given by 

Lfii/ = r,^ 11 ' — To 9fiv 
The notation 

(b,x) = b li x fl , ~g\b,b) = r v b»b v 

has been used. Associated to the pair (a: M (r), 6 M (r)), there is a preferred family of conformal 
factors, 6, which are developed along the congruence using 

iPV^e = 6(6, i), 61,5=6*. (12) 

Consistent with equation |l]), define g^ v = Q 2 g^, and denote the corresponding Levi-Civita 
connection by V. In this conformal gauge the 1-form along the conformal geodesies is given by 

/ = b-Q^dQ 

and satisfies (/, x) = 0. In the following, we will work with the unphysical metric g^ v and the 
1-form 

= e/ M + v M e. (13) 
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The 1-forms and /p are used to define a Weyl connection, V, via 



MP MP V MP 1 " p 

Vp + (^/ P + ^ P /p - Spp^/a) • (14) 



In the gauge corresponding to the connection V the equations (11a I and (lib I take the simpler 
form 

x^VuX" = 0, L A pi A = 0, 
where L pv is the Schouten tensor of the Weyl connection V. It is given by 

Note that generally one will have that L\u,v] ^ 0. 

We transport an orthogonal frame e^, k = 0,1,2,3, such that e = x, along the congruence 
of conformal geodesies according to 

± M Vpe fc = 0, (15) 

which implies that the frame is automatically g-orthonormal. If g is a vacuum metric, then the 
conformal factor, 0, can be written down explicitly by 

e^e, + e,(T-n) + ^e,(T-T,) 2 , (ie) 

where 0*, 0* and 0* are functions which are constant along a given conformal geodesic. They 
are subject to the constraint 

20,0* = g* (d», d,). (17) 
Furthermore, along each conformal geodesic 

d a = (d,e a ), 

are constant for a — 1,2,3 — see e.g. [17] . 

Analogous to the case of metric geodesic congruences there exist deviation vector fields for 
congruences of conformal geodesies [T7], referred to as conformal Jacobi fields and denote them 
by rf . Let r)k = ?7pe£. These fields satisfy 

V±?7fc = g(V v x,e k ), (18) 
We say the congruence has a conjugate point if r\ a (a = 1, 2, 3) vanish at that point. 

3.2 Conformal Gaussian gauges and space spinors 



The coordinates x^ described in section 2.2 will be extended off a fiduciary hypersurface by 
dragging them along a congruence of conformal geodesies. Hence, setting x° = r one obtains a 
coordinatisation of a slab of the fiduciary hypersurface. In these coordinates one has that 

e = x = d T . (19) 

Let L/ fe and Lij denote, respectively, the components of T " p and Lp„ with respect to the frame 
efc propagated according to the rule (15 1. One then has that 

fy' fe = 0, L ok = 0, (20) 

— see e.g. [T31 HI]. This choice of coordinates, frame field and conformal will be referred to 
as a conformal Gaussian gauge system based on the fiduciary hypersurface. The appearance of 
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conjugate points in the congruence of conformal geodesies implies a breakdown of the conformal 
Gaussian gauge. 

Let r M = \/ T lx ti , and let t aa be its spinorial counterpart. In terms of spinors, the gauge 
conditions (19 1 and (20 1 can be written as — see [IS] for details of the definition of Qaa'BB' — 



r AA 'e AA , 



V2d Tl t aa T 



AA' ' 'c — 0> ®AA'BB' = 0. 



The frame vectors eAA' will be expressed with respect to the vectors {co, ci, C2, C3} where cq is 
the normal to the surfaces of constant r and c\, ci and C3 are given by (8al-(8cl. We will write 



e A A' = e A A' c s: s = 0,1,2,3. 



In order to obtain a suitable set of evolution equations we consider a space-spinor formalism 
— see e.g. |30l 1181 127] — as it allows us to obtain evolution equations involving only unprimed 
spinors. Introduce a spin dyad {o, 1} such that 



T AA ' = A d A ' + i A l A ' . 

We decompose the connection as follows 

V A a< = \t aa ,X7 - t c a ,X7 ac , (21a) 
V AB =T {A B 'V B)B ,, V = r cc "Vcc<. (21b) 

Note that V AB is not the Levi-Civita connection of the surfaces of constant r, but the so-called 
Sen connection. The frame fields e s AA , are similarly decomposed using 

<?aa> = g 4' r cc ' t aa , - t c a , e AC , (22a) 
e AB =r (A B 'e s B)B ,. (22b) 

The fields 

e AB are associated to spatial vectors, and hence, they satisfy the reality conditions 

e AB — —t a a t b b e A , B ,. (23) 

Contracting with t aa and symmetrising, one decomposes the remaining connection and curva- 
ture spinors as follows. The space spinor Q AB cd — &ab(CD) is defined by 

QABCD = QaA'CC'T A B T C D — Q(AB)(CD) + ^AB^G '(CD)- 

One observes that Taa'SC = ^aa'BC+^abIca' and t aa f AA < = 0. Define T ABC d = t b b T ab > C d, 
which in turn, will be decomposed as 

FaBCD = (£,ABCD - X(AB)Cd) ~ ^ABfcD- 

The spinors in the latter equation possess the following symmetries 

^ABCD — ^AB(CD), XABCD — XAB(CD)i £,ABCD = £,(AB)(CD)- 

3.3 The conformal Einstein evolution equations 

In the sequel we will make use of propagation equations which are derived from the extended 
conformal field equations of [T3] — see also [THl HE] • These equations are a generalisation of the 
original conformal equations which allow the use of Weyl connections, as defined in (14). Weyl 
connections make it possible to consider more general gauges in the derivation of systems of 
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propagation equations out of the conformal field equations. In particular, it makes it possible to 
make use of the conformal Gaussian gauge systems discussed in section |3.2| 

Using the space spinor formalism and suitable contractions of the extended conformal field 
equations with t bb , one finds that the conformal Gaussian gauge system implies the following 
propagation equations for the unknowns e s AB , £abcd, Iab, X{AB)cd, ©(ab)cd and & G G CD : 



dre° AB = -X [AB) EFe %F ~ Jab, (24a) 

9 r e AB = -X {AB ) EFe EF, r = 1, 2, 3 (24b) 

d T £,ABCD — ~X( A b) EF £EFCD + —^(^ACX(BD)EF + ^BDX(AC)Ef) f EF 

~^X(AB)(C E fD)E ~ \{tAC®F F BD + e BD<d F F Ac) ~ '&VABCD, (24c) 

OtIaB = ~X {A B) EF fEF + -^0 F F ABl (24d) 

d T X(AB)CD =■ -X( AB ) EF Xefcd - Q(cd)ab + Qvabcd, (24e) 

d T Q(AB)CD = -X(CD) EF ®{AB)EF ~ d T Qrj A BCD + iV2d E \ a Hb)cde, (24f) 



r) Pi G - v EF 9i G 
VtVq ab — X(AB) U G EF 



V2d EF 



VABEF, 



(24g) 



where t\ A bcd and [iabcd denote, respectively, the electric and magnetic parts of 4>abcd and 
d A B is the spinor representation of d AA < — for further details see e.g. [TH [TH1 HZ] . Note that the 
above equations involve an unspecified conformal factor and the space spinor representation of 
the 1-form d^ (d A s) introduced in (13 1. For we will use (12 1 and (16 1. 

The evolution equations for the spinor 4>abcd are derived from a space-spinor decomposition 
of the Bianchi identity 

V AA '4>abcd = 0. (25) 
The resulting propagation equations ar e no t unique as one can always add to them a multiple 
of the constraint equations implied by (251, ~SI AB 4>abcd = 0. In [13j [14] the following Bianchi 
propagations equations 



V2d T (j)ABCD - 2V (£) 4>ABC)F = 



(26) 



are considered — to the so called standard system. Note that equation (26 1 is expressed with 
respect to the Levi-Civita connection V and not the Weyl connection V. 

In order to keep track of the possible appearance of conjugate points in the solutions to the 
propagation equations, we append to (24a)-(24gl and (26 1 an evolution equation for the Jacobi 
field. The conformal Jacobi field rf L has a spinorial counterpart T)aa' which can be split as 

VAA> = ^V T AA' - T B A ,rj AB , Tj = T) AA >T AA ' , J] AB EE T {A B T) B)B , . 



Conjugate points in the congruence of conformal geodesies arise if r/ AB = 0. Equation ( 18 ) implies 
that the fields 77, rj AB satisfy the propagation equations 



V2d rV = f AB v AB 



V2d- 



rVAB = XCD{AB)V 



CD 



(27a) 
(27b) 



Besides the above propagation equations ( 24a |-( 24g I , and ( 26 1 one also obtains a set of equa- 



tions referred to as the extended conformal constraint equations involving the operators cab and 
V ab (respectively V^s)- It should be noticed that despite their name, the constraint equations, 
are not intrinsic to the hypersurfaces of constant r — except, possibly, for an initial fiduciary 
slice — as generically, they contain d T derivatives. The crucial observation for our purposes con- 
cerns their propagation properties as given by the following lemma proved in |13j . Let D + (S) 
denote the future domain of dependence of S — see e.g. [22j [32] for definitions. 
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Lemma 1. Assume that the constraint equations are satisfied on an initial hyper surface S. If the 
conformal propagation equations are satisfied on D + (S), then the conformal constraint equations 
are also satisfied on D + {S). 



3.4 Structural properties of the evolution equations 



We discuss now some general structural properties of the equations (24al-(24g|, (26 1, ( 27a )-( 27b I 
which will be used systematically in the sequel. Let 



Vi = ( P(ABCD) i , 

where the subscript (ABCD)i indicates that after symmetrisation i indices are set to 1. Introduce 
the notation 

v = ( e 3lB , ^ ABCD , &ABCD , T)i VAb) , = (00,01,02,03,04), 

where it is understood that v contains only the independent components of the respective spinor 
— which are obtained by writing linear combinations of irreducible spinors. In terms of v and (f>, 



the propagation equations (24al-(24gl and (27a)-(27bl can be written as: 



d T v = Kv + Q(v, v) + Lcj), 



(28) 



where K and Q denote , respectively, a linear constant matrix-valued function and a bilinear 
vector-valued function both with constant entries and L is a linear matrix- valued function with 



coefficients depending on the coordinates. Similarly equation (26) can be written in the form 

V2Ed T cj> + A AB e\ B d- s 4> = B{T ABCD )<t>, 



(29) 



where E denotes the 5x5 identity matrix and A AB e AB , s — 0, ... ,3, are 5x5 matrices depending 
on the coordinates, while B(T abcd) denotes a constant matrix-valued linear function of the 
connection coefficients T abcd- For later reference it is noted that 
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and that 
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with r — 1, 2, 3. From the reality condition (23 1 one has that 



p° 



-p° 



'01 J 



'01 



'01' 



so that in particular and are real. The Hermitian matrices \[2E + A AB e AB , A AB e r AB 
imply real symmetric matrices for the propagation equations obtained from splitting the system 



( 28 1- ( 29 1 into real and imaginary parts of the propagation equations. These matrices are explicitly 



given in [27 
Let u = 



for some N E N. In 



(Re(f ), Im(u), Re(0), lm(0)). The unknown u takes values in 
this article u will be defined over S 3 or subsets thereof. From equations ( 28 1 and ( 29 1 it follows 



that u satisfies a system of quasilinear partial differential equations for u of the form 



3 



A u (u) + A r {u) ■ cp(u) + B{t, rr 4 , u) 
=i 



(31) 



10 



with Cf{u) denoting the vector fields (8a 
the matrix valued functions A s (z), s 



(8c I acting on the unknown u. 



Given any z E 
0,1,2,3 have entries which are polynomial in z. These 



polynomials are at most of degree one and have constant coefficients. The matrices are symmetric 
t (A s {z)) = A s (z), z E R N . The matrix valued function B = B{r,x A ,z) with {t,x a ,z) € 
Rx§ 3 x has entries which are polynomials in z (of at most degree 1) with coefficients which 
are analytic functions on K x § 3 . 



4 Purely radiative reference spacetimes 

In what follows we will construct purely radiative vacuum spacetimes using the ideas of [10 . These 
spacetimes are required to allow a conformal compactification such that the resulting conformal 
metric is smooth at null infinity including the point past or future timelike infinity. The purely 
radiative spacetimes in [TD] have been obtained from of static, asymptotically flat solutions to the 
Einstein field equations. Before discussing the results in [TU] which are crucial for our analysis, 
we introduce some definitions. 



4.1 Asymptotically Euclidean hypersurfaces 

Asymptotically flat spacetimes, {M,g^ u ), can be obtained as the development of asymptotically 
Euclidean initial data sets. The initial data (S ,h a p, K a p) will be said to be asymptotically 
Euclidean if there is a compact subset of S whose complement for some positive number ro is 
diffcomorphic to {y a E M 3 \ \y\ > r }, where \y\ 2 = (y 1 ) 2 + (y 2 ) 2 + (y 3 ) 2 - In terms of the 
coordinates y a introduced by this identification a standard asymptotic flatness requirement is 

h aj3 = - M + S a p + O (j^j ) as 1 2/1 -* °°> 

with to a constant — the ADM mass of S. An initial data set (5, h a p, K a p) is said to be time 
symmetric if K a p — 0. It follows form the Einstein field equations that its development will have 
a time reflexion symmetry with respect to the surface S. 

As we shall be working with conformally compactified spacetimes, it will convenient to consider 
a compactified version of the initial hypersurface S . To this end, it will be assumed that there 
is a 3-dimensional, orientable, smooth compact manifold (S, h), a point i E S, a diffcomorphism 
$ : S \ {i} — > S and a function tt E C 2 (S) n C°°(5) with the properties 

Q(i) = Q, D a Q(i)=0, D a D fl(i) = -2h a p(i), (32a) 
fl > on S\ {i}, (32b) 
h a p^n 2 ^^~h a p. (32c) 

The last condition shall be, sloppily, written as h a p — fl 2 h a p — that is, S \ {i} will be identified 
with S. Under these assumptions (S,h a p, K a p) will be said to be asymptotically Euclidean and 
regular. Suitable, punctured neighbourhoods of the point i will be mapped into the asymptotic 
end of S. 

In the sequel, we will consider initial data sets for the Einstein field equations which satisfy 
the boundary conditions (32a)-(32c) but have vanishing mass (to = 0). It follows from the mass 



positivity theorem that if to = then either (S, hup) is data for the Minkowski spacetime, or the 
initial data is singular somewhere in S — see e.g. [551 HH]- O ur discussion of initial data sets 
with vanishing mass will be concentrated on a small neighbourhood B a (i) of the point i where 
the relevant objects are suitably smooth. 



4.2 Time symmetric initial data sets with vanishing mass 

The construction of purely radiative spacetimes of [10] can be summarised as follows. Let V = 
(^^) 1//2 b e t ne norm of the timelike Killing vector, of some static, asymptotically flat, 
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vacuum spacetime. Assume that the spacetime has non-vanishing mass. In a static spacetime, 
there are coordinates (t,x a ) for which the spacetime metric takes the form 

-static = V 2 dt 2 + y^j^dx^x 13 , 

where h^ ttc = V~ 2 Zf a p denotes the metric of the hypersurfaces of constant t, whereas is 
the metric of the quotient manifold in the case of static spacetimes the quotient manifold 
and any surface of constant t can be identified. Denote by S one of these (time symmetric) 
hypersurfaces of constant t. The metric h^ tlc satisfies on S the time symmetric Hamiltonian 
constraint r [/j sta * IC ] = o, whereas 7 Q ^ and V satisfy the vacuum static Einstein field equations. 
Static spacetimes are usually analysed in terms of the quotient metric 7 Q| a. The asymptotic 
behaviour of the quotient metric 7^ is best understood in terms the behaviour of a conformally 
rescaled quotient metric 7 Q( a = (f2 statIC ) 2 7 a/ 3 in a neighbourhood, B a (i) C S, of the point at 
infinity i. As before we have 5 = 5U {i}. Using the so-called conformal static equations Beig & 
Simon [5] have shown that there is a choice of £l statlc and coordinates for which £l statlc and j a p 



are analytic in B a (i) and satisfy the boundary conditions (32a!-(32cl. The crucial observation in 
[TU] is that the conformally rescaled metric 



h aP = l -(\ + vfv- 2 ht¥ l \ 

also satisfies the time symmetric Hamiltonian constraint. The metric h a p is asymptotically 
Euclidean with vanishing mass. Let 

n = i(i + v) 2 v- 1 n static , h a0 = (n staUc ) 2 h Ql3 = n 2 h s Js tie . 

v2 



The pair (h a p,£l) satisfies the time symmetric conformal Hamiltonian constraint, equation (6a I 
with £ = . Bo t h and h a /3 are analytic in B a {i) and, furthermore, they satisfy the boundary 
conditions (32al-(32cl. The pair (h a R, fi) can be used to construct initial data for the conformal 
Einstein fie. 



d equations of section |3.3| — this data consists essentially of the value of the rescaled 
Weyl tensor and the Schouten tensor on the initial hypersurface. This initial data is obtained out 
of algebraic manipulations of the conformal constraint equations. Although h a p and f2 are given 
in a gauge which makes them analytic around i, the data for the conformal field equations need 
not be. It is a non-trivial result that the data constructed out of this procedure is indeed analytic 
around i. The time evolution of this data yields — using, say, the Cauchy-Kowalevska theorem — 
in a 4-dimensional neighbourhood of i an analytic spacetime metric g M „ and an analytic conformal 
factor O. 

In the sequel this local existence problem will be solved by a different method. Under these 
circumstances, the physical metric, g^ = 0~ 2 <? At! y, obtained on the timelike future I + (i) of i in 
this evolution is a solution of the Einstein vacuum field equations with an analytic structure at past 
null infinity, for which the point i represents a regular past timelike infinity i~ . These solutions 
of the Einstein field equations are radiative spacetimes in the sense discussed in the introduction. 
Note that due to the time reflexion symmetry, the same holds for the pair (/~(i), g^ u ) where i now 
represents the regular future timelike infinity i + . On the complement of J + {i) U J {i) — these 
are the points in the development which are spacelike related to i — one also has a solution of 
Einstein's field equations in a neighbourhood of i, for which the point i now represents the point 
spatial infinity i°. Definitions of the causal sets J* and can be found in [2"2l |3"2"] . 

The results of |10j relevant for our analysis are presented in the following theorem. 

Theorem 1. Given an analytic solution to the conformal static equations on B a {i) C S, there 
exists a solution, (h a p,Q), to the time symmetric conformal Hamiltonian constraint on B a (i) 
satisfying 

(i) it is asymptotically Euclidean and regular in a neighbourhood B a (i), 
(ii) it has vanishing ADM mass, 
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(Hi) the spinorial fields <Pabcd and Qabcd constructed out of (h a p,fl) are analytic in B a (i). 

It was long conjectured that given a series of multipole moments satisfying an appropriate 
convergence condition, there exists a static solution of the Einstein equations with precisely those 
moments. This conjecture was proved in full generality in [19] — see also |TJ [23] for partial results 
and an alternative approach. These results show that the metric, once written in a certain gauge, 
is uniquely determined by the multipole moments. Thus there is an infinite family of static 
metrics analytic in a neighbourhood of i. This family is parametrised by the multipole moments. 
As a consequence of Theorem [l] it follows that we have an infinite family radiative spacetimes 
constructed in the above mentioned way that will serve as reference spacetimes for our stability 
analysis. 

4.3 Construction of radiative reference spacetimes using conformal Gaus- 
sian gauges 

In this section we discuss how conformal Gaussian gauge systems can be used to construct radia- 
tive spacetimes out of the Cauchy initial data provided by Theorem [l] The original construction 
discussed in [TO] makes use of a propagation system derived from the "original" conformal field 
equations of [5JIH1IE]. As mentioned before, due to the analyticity of the setting, existence of a 
unique analytic solution in a spacetime neighbourhood of the point i follows from the Cauchy- 
Kowalewska theorem — see e.g. [4]. 

Here we follow a different approach. Our intention is to use the development of the data given 
by Theorem [l] as a reference solution from which a stability result will be derived. As seen in 
|27j . the discussion of the properties of the conformal boundary of radiative spacetimes simplifies 
with the use of a conformal Gaussian system. Thus, it is natural to discuss the structure and 
existence of the reference solution using the same type of gauges. In our approach, existence of a 
solution of the propagation equations implied by the extended conformal field equations and the 
conformal Gaussian gauge follows from a variation of a theorem of Kato — see [9| l26|, 127] . 



4.3.1 Setting up conformal Gaussian coordinates 



Let (h a p, ft) o n S be one of the solutions to the time symmetric conformal Hamiltonian constraint, 
equation (6a I, given by Theorem [l] It is assumed that £ = 0. In what follows all quantities 



derived from this initial set and from its time development will be distinguished with an overbar. 
Although it is customary to make use of a conformal factor ft which is non-negative, it will be 
convenient to consider a non-positive conformal factor on S — the reason for this will be evident 
in section [4. 3. 3| Accordingly, if one takes the conformal factor of Theorem [I] with an extra factor 
of —1, so that it is negative at least in a neighbourhood away from i. One has that ft satisfies 
the appropriate modification of the boundary conditions ( 32a l-( 32b I . Namely, 



ft<0onS\i, ft(i)=0, D A O,(i) = 0, D A D B Q(i) = 2h A B («)• 



(33) 



On a suitably small neighbourhood, B a (i) C S, we use the coordinates x A , A = 1,2,3 centered 
at i and consider the following initial data for a congruence of conformal geodesies: 



0, & 



e* = ft, e* = o, 4 = e*&* = (de)*. 



(34) 



The coordinates x A are extended off S by dragging along the congruence of conformal geodesies 



3.2 It follows that d* = at 



to obtain conformal Gaussian coordinates as discussed in section 
i and using d a (f) — d a (0), we get d a = D a (fl). The constraint ( IT I applied to our situation 



gives the initial data for O on B a \i, which must hence be positive there. The initial value at i 
is obtained taking the corresponding limit. Making use of normal coordinates centred in i, one 
finds that since ft = 0(\x\ 2 ), we have /i»(Dft,Dft) = C(|a;| 2 ) and thus § = 0(1). Hence §» > 



everywhere on S. It follows from equation ( 16 1 that along each conformal geodesic the conformal 
factor is given by 



0(f) = ft 



1 



0* 



= ft 1 - 



(35) 
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where 

uj = \ l is—, u)(i) = 0. 

V ©* 

Define now the conformal boundary, J? , in a natural way as the locus of points in the development 
of the data on S for which = and d0 ^ 0. It is easy to see that a conformal geodesic with 
data given by (34 1 passes through J? whenever f = ±uj. Note that every conformal geodesic, 
except for the one passing through i, intersects J 2 " twice, once to the past of S and once to the 
future. 

Remark. We shall always work on a neighbourhood B a (i) C S which is assumed to be 
sufficiently small to ensure h a P D a ClDpCl ^ in a neighbourhood of i. 

4.3.2 Existence and conformal structure of the reference spacetime 



We discuss now the existence of solutions to the propagation system (31 1 with data prescribed 
on B a (i) C S. The initial data for this initial value problem is obtained through the conformal 
constraint equations on S 0[1O] and satisfies on B a (i) C S 

e° AB = 0, (36a) 

e r AB = ®ab i (36b) 

Jab = 0, (36c) 

Cabcd = €abcd[7i], (36d) 

X{ab)cd = 0, (36e) 

Qabcd — —-^D( AB D CD )£l + —tlABCD, (36f) 
4>ABCD = ^D(AbDcD)^ + ^SABCD, (36g) 

r) = 0, (36h) 

VAB = &AB J ( 36i ) 

where Habcd = —£A(C e D)B and &ab are ^ ne spatial Infeld-van der Waerden symbols. The 
spinorial fields £,ABCD[h] and sabcd are the connection coefficients and the trace-free Ricci 
curvature of the 3- metric h a g . 



Since B a (i) is assumed to be a simply connected neighbourhood of i, we can find a diffeomor- 
ii 

for 



phism <f>B onto McS 3 . We extend the initial data to § 3 in the following way. The initial values 



e^B' fAB, X(AB)CD) V, VAB 

are constant and are extended with their respective values to all of S 3 . For 

6*, d*, Cabcd, Qabcd, §abcd (37) 



we use the linear extension operator E given in section 2.3 



to obtain fields defined over § 3 . The 



extended fields will be denoted on S 3 in the same way as in ([37 1. 

The vectorial unknown 

u(t,x) = fRe(u)(r, x), Im(u)(r, x), Re(^)(r, x), lm(0)(r, x)j 

is viewed as a function of r with values in the Sobolev space H S (§ 3 ,R N ). Let u (x) = u(0,x) 
denote the extended initial data on S 3 as described in the previous paragraphs. For 5 G R, m € N, 
m>2, define set 

Df = {w £ H m ($ 3 ,R N ) | (z,A°(w)z) > S(z,z),Vz ER N } , (38) 
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Figure 2: A conformal representation of the reference solution obtained as the development of 
Euclidean initial data with vanishing mass on S. The shaded region represents the physical 
spacetime. 



whe re A°(w) denotes the matrix valued function defined by the symmetric hyperbolic system 
(31 1, and (•, •) is the standard scalar product on M. N . It can be verified that 

A°(u ) = diag (l, . . . , 1, ±=, V2, V2,V2,^=, ±=, y/2, y/2, sfi, -^J . 

and hence u £ D™ for < 8 < l/y/2. 

One has the following local existence result. 

Theorem 2. Suppose m > 4. There exists a T > 0, a subset Df C D™ and a unique solution 
u(t) of equation (31) defined on [—To, To] with initial data uq on S ~ § 3 and such that 

u S C(-T 0) T ; DfJ n C 1 (0, T; H fr ^ 1 (S 3 , R N )). 

Furthermore, u £ C°°([-f ,f ] x § 3 ) and T can be chosen such that u(t) has non-vanishing 
J acobi fields t]ab, so that the solution is free of conjugate points in [—To, To]. 

Proof. The first part of the theorem follows from the generalisation given in [3] — see also [57]— 
of Kato's existence and stability result for quasilinear symmetric hyperbolic systems [26j . For the 
Jacobi fields we observe that r/^g satisfies an ordinary differential equation along the curves of 
the congruence. Hence there is a minimum interval for which the fields do not vanish. The lack 
of conjugate points then follows from our earlier discussion. □ 

Remark 1. It follows that on (—To, To) the solution, as given by the theorem, is of class 

H m ((-T Q ,T ) x § 3 ) c C m - 2 ((-T ,fo) x S 3 ). 



Remark 2. Let M = [—To, T ] x§ 3 and denote the metric on M. by g^ v . The spacetime (M, g^) 
will be used as the reference spacetime for our stability analysis later on. Let D(B a (i)) denote 
the domain of dependence of B a (i) — see e.g. [251 122 f° r definitions. The value of To can be 
chosen small enough so that 

M = Mn J-(i)cD(B a (i)). (39) 

The finite reduction of To is not a problem in our scheme as any To > corresponds to an 
infinite time interval in the physical description. As it will be seen in the sequel, it may be 
necessary to further reduce To by a finite (non-negative) amount. 

The spacetime (M \ J?, g^ v ) is conformally related to a vacuum spacetime, (M \ J?, 6 _2 <?^i/), 
with vanishing cosmological constant. The spacetime (A4\ J? , Q~ 2 g^) is a radiative spacetime as 
the conformal metric g M „ and associated fields extend smoothly through the conformal boundary 
J 1 and furthermore, there exists a point i + = (0, i) £ {0} x S C At satisfying (by construction) 
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the conditions to be a future timelike infinity — cfr. the discussion in the introduction. As a 
result of Lemma [T] it follows that the conformal constraints are only propagated in D(B a (i)). 
Consequently, the region [—To, T ] x S 3 \D(B a (i)) is of no physical interest as it does not describe 
a vacuum spacetime. It is nevertheless relevant to use the whole of (M^g^) for the correct 
formulation of the stability part of our work. In the following, in a slight abuse of language, we 
will refer to {M,g^ u ) as the purely radiative reference solution. 

4.3.3 Hyperboloidal data for the radiative reference spacetime 

The conformal affine parameter f defines, in a natural way, a foliation of the manifold M. Let 
S T denote the surfaces of constant f. For fixed f one has that S T is diffcomorphic to S 3 . Let 
r G (0,T ) and define 

S = {-f }x§ 3 , Z = {peS Q \Q = 0}. 

The set So intersects null infinity in a hyperboloidal way. More precisely, Z divides Sq into two 
regions: one where O < and one where O > 0. The latter contains the origin and will be 
denoted by H. Hence 

H = {peS o \e>0}. 

The latter definition justifies the choice of a negative conformal factor f2 made in the previous 
sections. We rewrite the conformal parameter f and the conformal factor of the reference solution 



given by equation (35), so that their initial surface is Sq. Define 

T = f + f , 0(t) = 0(t-t o ) 



so that r = on Sq and 



r 2 \ 1 



e(r) = n(i-^i- e,f T + -t*q, 



- 9 / — 9 ' 



The initial value of on Sq will be denoted by Cl. Recalling that t±(x a ) — ±lo{x a ) gives the 
location of null infinity, we can see that 

n>o M>|f |, 
n = o \Q\ = \f \, 

n<o <=> |w|<|f |, 

and hence Cl is the correct boundary defining function on So with Z = {p S So\Cl = 0} = &H.. 
One can briefly verify that the new conformal factor gives the correct location of null infinity. 
If Cl ^ we get that = whenever 



T ± UJ. 




Along the unique conformal geodesic 7 passing through the point N on S 3 at all times, we have 
f2 = and 

0= 1&*(t-t ) 2 . 

Thus, along 7 one has that = only when t — tq. This corresponds to the point i in the 
Cauchy hypersurface S. 

In what follows, let 

u(t, x) = u{t —Tq,x). 



16 



For consistency, we shall denote entries of u, as defined above, also with . The entries of u imply 



on So data for the symmetric hyperbolic system (31 ). We shall write Uo(x) = u(0, x). The data Uq 



are not truly hyperboloidal data. The reason for this is the following. The frame vector eo which 
on S is normal to this surface, ceases to be normal to the surfaces S f , of constant f if f ^ 0. This 
is a feature that differentiates conformal geodesies from standard geodesies. This phenomenon 
manifests itself in the fact that the components e AB of the space spinor decomposition of eo 
according to formulae (22a) and (22b), which are zero on S, become, in general non-zero off 5^ 



As a consequence, the spatial vectors of the frame e a , a = 1,2,3 which span the tangent bundle 
to S, TS, in general, pick up components off the tangent bundle TSr- Accordingly, the data ito, 
which is expressed in terms of the spatial frame e a (cab) is not intrinsic to So- We shall call this 
type of data tilted hyperboloidal data. 

In order to relate the tilted data uq to proper hyperboloidal data one proceeds as follows. 
Denote the normal vector field to So, the 3-metric induced by and the second fundamental 
form of Sq by n M , h a g and K a p respectively. Note that K a p and Xaji are two conceptually 
different quantities — the former is related to the hypersurface, the later to the conformal geodesic 
congruence. Choose a /i-orthonormal frame e , a = 1, 2, 3 and set e = n, so that e fe , k = 0, 1, 2, 3, 
forms a g-orthonormal frame on So- The frames {v,e a } and {n,e a } are related by a Lorentz 
transformation A* - = A l j(x), which will be regarded as a matrix valued function over § 3 , as 
follows 



A l ,A k 3 ^S\, (40a) 



j 

L,=A h j e k , (40b) 



-j 

The fields r.\, L^, dijki are the frame components of spacetime quantities and the Lorentz 

-i fc5 —ip —ijkl 



e,=A 7 fc e fe . (40c) 



transformation induces transformed components /., f/^, L^, daki m the canonical way. For 



example 



d — A m A n A p A 9 d 

yHjkl — yi i yi j yi fe yi l a mnpq- 



The SL(2,C)-spinorial counterparts of the above tensors can be obtained by suitable contractions 
with the constant Infeld-van der Waerden symbols. These, in turn, can be space-spinor split with 
respect to t m = v2e to obtain proper hyperboloidal data on Ti.. 



5 Radiative spacetimes from perturbed hyperboloidal data 

Given one of the radiative spacetimes discussed in section [4] it is natural to ask whether it is 
possible to establish for these spacetimes an analogue of the semiglobal stability results for the 
Minkowski spacetime discussed in [THIHT]. In this section we establish such a result. This result 
generalises the aforementioned stability results to the case where the reference solutions are not 
flat. 

5.1 On hyperboloidal data 

In order to discuss appropriately the notion of closeness of two hyperboloidal initial data sets, 
we make use of the notion of limit of spacetimes as discussed in e.g [21]. Assume that one has 
a smooth family of hyperboloidal initial data sets (Tt £ , h^p, K^p, f2 e , S e ), parametrised by e > 0, 
such that 

(H ,h a/3 ,K a p,n ,£ ) = (H,h a 0, K a p,&, £), 
where TL e are simply connected 3-dimensional manifolds with boundary &H e . The fields {h e a p, K^g, fl e 



satisfy the vacuum constraints, equations (6a) and (6b I, with 



il £ = 0, dn £ ^ on dH e , n £ > on the interior of H £ 



^^Note, however, that in the case of the conformal Minkowski spacetime in the form discussed in |27) . tangent 
vectors to the congruence of conformal geodesies are always normal to the surfaces of constant affine conformal 
parameter. 
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As the manifolds TL e are taken to be simply connected, there exists a family of diffeomorphisms 

<ff '. T~C e — ► IA C § 3 , <j) e {dH e ) = dU, 

with IA as in section [2. 2 1 In particular, 

$ £ = 4)- 1 o cj) £ : H £ -> W with = 0°, 

is a smooth family of diffeomorphisms between 7i e and H, such that <I> is the identity map. Note 
that the locus of points for which J7 e = coincides as sets on S 3 , by construction, for all e > 0. 

Let {e^}, a — 1, 2, 3 denote a smooth family of ft, e -orthonormal frames such that 

e° = e„, a = 1, 2, 3, 



with e a as given in the end of section 4.3.3 The Levi-Civita connection D £ of the metric h e a p, 
together with the frame (f a , imply the 3-dimensional Ricci tensor r[h e ] a p and spin connection 
coefficients 7 a h c , as well as their spinorial counterparts r ABCD and £ ABCD - From these quantities 
and the tensor K e a a one can determine the e a -components of the Schouten tensor and e -electric 
and e -magnetic parts of the rescaled Weyl tensor for the spacetime on H £ — see e.g. [H [§]• 
Using the frame e a and r M = \/2g , their spinorial counterparts can be readily calculated. This 
procedure produces proper hyperboloidal data — in the sense discussed in section |4.3.3[ We note 
that so far the data is only given on Tt E , respectively, (j) £ (Ti s ) = Wc§ 3 . We use the extension 



operator E given in section |2.3| to extend the proper hyperboloidal data to the whole of S 3 , 



bearing in mind that outside IA the data may violate the vacuum constraints ( 6a I and ( 6b ) 



In the following, the parameter e will be dropped from the expressions, in order to ease the 
presentation. 

5.1.1 Local existence for hyperboloidal data 

From studying the proper hyperboloidal data on the hypcrsurfacc Sq in the reference solution 
one can see that this data does not include all the quantities that are needed for the tilted 



hyperboloidal data used for (31 1. We are missing the data for e s AB , /ab, £abcd, Xabcd- 

In order to consider the construction of tilted hyperboloidal data for the conformal propagation 
equations we shall first develop a small part of the perturbed spacetime. For this we need a 
local existence result for the development of data constructed using an auxilary congruence of 
conformal geodesies which departs the initial hyperboloid orthogonally. The information thus 
obtained complements the proper hyperboloidal data constructed above. 

For the construction of the auxilary congruence we use the following initial setup. Let x = n, 
ek\s = e k , i ABCD \s = S,ABCD[h], X {AB)CD \s = K ABC d- As long as the data for f_ AB is smooth 
and satisfies (/, n) = 0, i.e. f AB — f t AB \ ^ can be chosen arbitrarily, since we are only interested 
in local cxistience in some small neighbourhood of the initial surface. In particular, we can set 
■Lab = a S ain ' 

We can now proceed to our local existence result for the unknown 

— = ^ab i &AB , Lab ' ^abcd ' *(ab)cd ' ^-Abcd > <t AB cD ' 1' ^ab^ 

The underline in the above spinors indicates that they are expressed with respect to the frame 
e a . The proposition below is essentially very similar to Theorem[2] This is not a stability result, 
but purely a local existence for the perturbed data. For this we define DJ 1 analogously to ( 38 ) . 



Proposition 1. Suppose m > 4. There exists a T > 0, a subset ]J T C D™ and a unique 
solution u(t) of equation (31) defined on [— T ,T ] with initial data u on S ~ § 3 and such that 

«eC(-T o ,T ;Dr o )nC 1 (0,r; ff™- 1 ^ 3 ,^)). 

Furthermore, uG C°°([— T , T ] x § 3 ) and u(r) has non-vanishing Jacobi fields r\ . , so that the 
solution is free of conjugate points in [—Tq,Tq\- 
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Figure 3: The auxilary congruence (— — — ) normal to H. and the tilted congruence ( — ) used to 
evolve the perturbed initial data. 

Remark: The local existence theorem clearly applies to (it, h a R, K a f3: and the resulting 
spacetime coincides locally with the reference spacetime (Ai^g^). One could then consider the 
derivation of a result similar to the stability of Minkowski in [57] . However it should be noted 
that there is no guarantee that this setup covers i + or J?+ of (Mjg^v) and its perturbations. 
Thus a semiglobal existence and stability result can not be guaranteed. 

5.1.2 Construction of tilted hyperboloidal data 

As a result of Proposition [T] we have a spacetime slab TV = [— T ,T ] x § 3 with metric g^ v which 
describes a vacuum spacetime on A4 = J\f f] D(H). For sufficiently small e, one can regard the 
3-manifold S as a hypersurface of a spacetime (TV, g^v)- On this hypersurface we now want to 



give tilted hyperboloidal data for (31 1. For this we will evolve a conformal geodesic congruence, 



that departs S in the same way as in the reference solution. Thus we can locally obtain all the 



components of the unknown u on S and use this as initial data uq for (31 1 



We know fix the initial data for the frame and the congruence. As before, let n denote the 
unit normal to S. Then {e^} = {n,e a } is an orthonormal tetrad on S. From e^- we construct a 
tilted orthonormal frame e^ on TJ\f\s via 



where A^- is the same matrix valued function associated to the Lorentz transformation via (40a 



(40c). It follows then that the timelike vector eo has the same projections on the orthonormal 
frame e_j as eo has on e,-. That is g(eo,e_j) — g(eo,e_j)- Note that eo is by construction tangent 
to a congruence of conformal geodesies. For discussions of the development of the perturbed 
initial data (H.,j a p, K a p,Sl,Y,) which go beyond local existence, it will be important to consider 
a congruence of conformal geodesies that departs the hypersurface TL in the same way as the 
congruence of conformal geodesies in the reference spacetime arrives at Ti.. For this reason we 
will use x = eo. 

As in the case of the reference spacetime (A4,g^), we want to cover the development of the 
initial data set (H, h a p,K a 0, fi, E) with a conformal geodesic coordinate system. The formulation 
of the evolution equations is dependent on a gauge for which (/, v) = 0. Taking this as well as 
d^ = O/jn + V^O into account, we have 



do = (d,e )=A I] + A a e a (tt), (41) 
d a = (d, e a ) = Qf a + A a °S + A a "e fc (SI) . (42) 

Assume that the perturbed hyperboloidal data (Ti., h a p, K a p, SI, S) satisfies the technical condition 

d a ~A a Z-A a b e b (n) = O(Sl), (43) 

close to &H, where d a denote the value of the components of the 1-form d^ for the reference 
solution and 0(SI) indicates that the quotient of the left hand side divided by SI goes to zero as 
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f2 goes to zero. Initial data for the congruence of conformal geodesies on which the conformal 
Gaussian system hinges will be prescribed such that: 



x = e , 
d„ = d„ 



f a = n- 1 (d a -A a °X-A a b e b (n) 



(44a) 
(44b) 

(44c) 



on TL. As a result of the technical condition (43 1 one has that that f a has a finite limit at &TL 
— and is thus well defined there. This particular choice of initial data for the congruence of 
conformal geodesies provides an easy analysis of the structure of the conformal boundary of the 
spacetime — as will be seen in |5.2| 

Initial values for 0*, G*, G* which characterise the canonical conformal factor G associated 



to the congruence of conformal geodesies defined by (44a)- (44c) are given on TL by 



6* = fi, e* = (d,e ), 2ftG* =g i (d,d)*, 



(45) 



with the corresponding limits for points on dTL. As usual, the subscript * indicates that the value 
of the functions are extended of the initial hyperboloid by requiring them to be constant along 
a given conformal geodesic with initial data given by (44al-(44c). It follows from the conformal 



Hamiltonian constraint, equation (Gal that d is null on dTL. Note that as a consequence of the 



conformal Hamiltonian constraint it follows that when fl = 



D a rtD a n = o, 



consistently with the choice of initial data for the conformal factor in (45 



Initial data for the Jacobi field associated to the congruence of conformal geodesies arising 

Vk = Vk = {v>h), onH. 



from (44a I- (44c I is set by 



Finally, one needs to construct tilted data for the rescaled Weyl (d^xp) and Schouten (L^) 
tensors. For this, one starts with the standard hyperboloidal data for these tensors implied 
by the conformal constraint equations. As a consequence of the local existence result given by 
Proposition [TJ one can regard TL rightfully as a hypersurface of a spacetime. Let dyjy and 
denote the components of the Weyl and Schouten tensors with respect to the frame {e k }. The 
tilted initial data for the Schouten and rescaled Weyl tensors on TL is given by 

Lij = A/ l A J m L nm , diju = A i m A :( "A /i , p Aj 9 d mnpg . 

The required spinorial data is then obtained by means of a space spinor decomposition of the 
spinorial counterpart of the above tensor. 



5.2 Structure of the conformal boundary 

In this subsection we analyse the structure of the conformal boundary of the development of 
hyperboloidal data which is close to hyperboloidal data of the reference radiative spacetime. 
Attention is focused, in particular, on the location of time like infinity i + . This investigation is 
possible as in the vacuum case the conformal factor can be calculated a priori from the initial 
data. 



Along the congruence of conformal geodesies with initial data on H given by (44a)- (44c I, the 
canonical conformal factor, G, takes the form 

0(r) = n + Q,r+ \&*t 2 : (46) 
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with 0* and 0* given by the conditions in (45 1. As it is customary, denote by J> the locus of 



points in the development of the data (hij,Xij, 0, S) for which = 0. For q S J" to qualify as a 
candidate to be the timelikc infinity of the development it needs to satisfy 

0(g) = 0, dO(q) = 0, Hess0(g) non-degenerate. (47) 

We recall that the 1-form d M = 9/^ + V M and that 200 — g^(d,d) hold. Thus, for smooth / M 
and O we observe that on J* , c? M is null and can only vanish there if d0 — 0. Therefore the only 
point on at which all d a — (d, e a ) can vanish, is at a candidate for i + . 



Define the function 



A = g\d*,d*) - Ol = -S ab d a d b , 



and observe that A = if and only if d a — 0. We recall here that d a , and hence A, are constant 
along conformal geodesies — see e.g. P^l [T7] . 

Now, by construction of the reference solution from the initial Cauchy data given on S and our 
choice of B a (i), there is a unique point in B a (i) C S such that A = 0, namely i itself. By our choice 
of coordinates, i lies at the origin of our spatial coordinate chart and is identified with the north 
pole of S 3 . Thus by Theorem bland remark 2 following it, the point (f,x A ) = (0,0,0,0) is the 
unique point in D(B a (i)) satisfying the conditions for i + . This corresponds to (r, x A ) = (fo, 0, 0, 0) 
in the 5o-adapted coordinates. Furthermore the points (r, 0,0,0) are by construction the only 
points in D(B a (i)) where A = for the reference solution. Given our choice for the initial data 



of dfj, in (44b I for the perturbed spacetime, it follows immediately that the only point at which A 
vanishes on TL is the origin. Therefore the only candidate for i + is the point, assuming it exists, 
where the conformal geodesic 7(7-) — (r, 0,0,0) intersects . Note that for general initial data 
it is possible for the development to contain a singularity before reaching i + — as in the case of 
black hole spacetimes. 



The first condition in ( 47 1 implies 



(48) 

while the second one can be rewritten as 

dtt + dO*r + ^dO*r 2 + (0. + 6*r)dr = 0. 
Together with A = 0, both conditions imply that 

r = T+s -^ = -4L. (49) 

For sufficiently small perturbations, (9* = (d, v) = A °S + A a e a (f2) does not vanish and r + is 
finite. In fact later on the initial data will be restricted such that t + lies in the interval (0,2Tb). 
Note that for the reference solution we have r + = T and q = i 6 S . By construction we have 
that Hess0(<;) is non-degenerate, and hence the i + -candidate in the reference spacetimes can be 
regarded, rightfully, as future timelikc infinity. 



5.2.1 Behaviour of the Hessian of O 



The point q = (r + , 0,0,0) satisfies the first two conditions of (47 1. We now verify that the 
conditions discussed in the previous paragraphs imply that the Hessian of is non-degenerate at 
q. For this we make use of the following expression: 



V„V y = V.V^ + VuV(,0,t 



2V (/1 0*V^T + 2rVf u 0*V^T 



-(©* + eUrJV^Vi/r + 0*V AI rV, y r. 
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Using the conformal Gaussian coordinates (r, ar 4 ) for our calculations and a = (b,v)* = 0*/f2 as 
well as ( 49 1 we find that 



V V 6| l+ = W, 



2n. 



(50a) 
(50b) 
(50c) 



where the right hand side is given in terms of initial data and spatial derivatives. Observe that 
V Ada are identical to the ones in the reference solution and thus the only change arises from the 
perturbation of a. For data 0*,f2 close enough to the reference solution, in the sense that at 
the spatial origin a and V_^a are sufficiently close to their counterparts in the reference solution, 
the Hessian will not be degenerate and thus q can be identified as the timelikc infinity i + of the 
perturbed spacetime. Note that one may rescale the initial data such that X remains fixed, while 
il = Ct. This greatly simplifies the analysis of the effect of perturbations of a on the Hessian. 



5.3 The stability result for purely radiative spacetimes 

We are now in the position of presenting our main result — a stability result for purely radiative 
spacetimes arising from hyperbolidal data. The proof of this result follows the ideas of [101 [27] . 



Let uq and uq denote, respectively, the tilted initial data for the radiative reference spacetime 
and for a perturbation thereof, as discussed in the previous section. Using the diffeomorphism <p, 
we can think of Uq and Uq as a vector- valued function over S 3 . On the other hand, using </>, u 
can only be regarded as a vector valued function over Wc§ 3 . On U, let 



W 



U - U 

Using the linear extension operator given in section 

Uq = Ew, Uq = 



2.3 



with Ew satisfying (10b I, we set on S 
m + Uq. 

Consistently with the previous discussion of the initial data, one writes 

u = u + u, (51) 

where we recall that both u and u are regarded as vector valued functions over K x § 3 . In 
particular, the quantities in 

u = i&AB j fAB 7 £,ABCDi X(AB)CD > ® ABCD i 4>ABCD j VAb)i 

describe the (non-linear) perturbations of the reference radiative spacetime. Substituting the 
Ansatz (51) into the evolution system (31 1 one obtains a symmetric hyperbolic system for u of 
the form 



3 

A°(u + it) ■ d T u + ^2 Af (u + ") ' c p(") + ^( T ' x " 4 ' u,u)-u = 0, 

r=l 



(52) 



with B(t, a; -4 , u, it) a matrix valued function with entries which are polynomials in u of at most 
degree one and coefficients which are smooth functions on K x S 3 . As discussed earlier, the matrix 
valued functions A s (z) are symmetric, have entries which are polynomial in z of at most degree 
one and have constant coefficients. In contrast to the situation discussed in |27j . the matrix 
A°(uq) is not diagonal. Nevertheless, Theorem [2] guarantees that there exists a neighbourhood 
Df a C D™ such that uq £ Df Q . Hence, 



(z, A°(u )z) > 5(z,z), Vz e 



niV 



The previous discussion allows us to formulate an existence and stability result for the solutions 



to equation ( 52 1 
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Figure 4: The perturbed radiative spacetime and its conformal boundary. Note that this figure 
is not a conformal diagram. 



Theorem 3. Suppose m > 4. Let Uq — uq + Uq be hyperboloidal initial data satisfying the 
technical condition (4-3). Given Tin(fo, 2fo), there exists e > such that: 

(i) For |J tt nJJm < e t here exist a unique solution u — u+u to the conformal propagation equations 
(24-a)-(24g) and (26) with minimal existence interval t S [0,T] and u G C m ~ 2 ([0,T] x § 3 ). 

(ii) The associated congruence of conformal geodesies contains no conjugate points in [0, T]. 

(Hi) At the origin one has t + = — £1/0* G [0,T] . 



(iv) The Hessian as given by (50a)-(50c) is nondegenerate at (r+, 0,0,0). 

The solution u = u + u on T> + {S) implies a C m ~ 2 solution (Ai,g) to the vacuum Einstein field 
equations with vanishing cosmological constant, where g^ = Q~ 2 g^ v with © given by (46). 



Proof. By hypothesis, we know that uq G Df c D™. We thus satisfy the conditions of 
the variation of Kato's theorem given in [101157]. Parts (i) and (ii) follow directly from Kato's 
theorem. Parts (iii) and (iv) follow from the discussion in section 5.2 Due to Lemma [l] on the 
propagation of the extended conformal constraints, one knows that the existence of a solution of 



the propagation system (24al-(24g) and (26 1 implies a solution to the full extended conformal 



Einstein field equations which in turns implies a solution to the vacuum Einstein field equations 
with vanishing cosmological constant. □ 



One has the following corollary. 

Corollary 1. Suppose the condition of Theorem^are satisfied, then a solution (A4,g), as given 
above, has a conformal boundary given by the set = 0. The conformal boundary consists of 
the set , which represents future null infinity, and the point i + given by (t+, 0,0,0), which 
represents timelike infinity. 

Hence the conformal boundary of radiative spacetimes is shown to be stable, subject to the 
conditions of Theorem [3l 



Remark. Note that the evolution of uq implies a unique solution u and a spacetime, which is 
locally isometric with (TV", g^) as obtained in Proposition [T] 



6 Conclusions 

The conformal Einstein field equations have been used to formulate an hyperboloidal initial value 
problem by means of which one can address the nonlinear stability of a class of purely radiative 
spacetimes. Our analysis is particularly concerned with the structure of null infinity and the 
location of future timelike infinity, i + . The use of a conformal Gaussian gauge system coordinates 
allows to calculate the location of the conformal boundary directly from the initial data, as well 
as analyse its structure. The results are a generalisation of previous results on Minkowski-like 
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spacetimes [27 . In particular the conditions on the location of i + are given in a more general 
form, while the choice of initial data for the perturbed conformal geodesic congruence facilitates 
the stability analysis of the conformal boundary in comparison to [27 . 

The reference purely radiative spacetimes are constructed from static initial data using pro- 
cedure introduced by Friedrich [10]. Recent results on the necessary and sufficient convergence 
conditions required on a series of multipolar data, one can conclude that there is an infinite family 
of reference purely radiative spacetimes. 

We observe that we were able to prove the stability of radiative spacetimes and study their 
conformal boundary without requiring any explicit knowledge of the reference spacetime. All our 
analysis has been carried out in the abstract only relying on general properties of these spacetimes. 
This raises the hope that similar results may be possible for other classes of spacetimes. 
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